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I— I A deterministic theory describing the behavior of an ultrathin lubricant film between two 

Q atomically-smooth solid surfaces is proposed. For the description of lubricant state the param- 

C^ eter of excess volume arising due to chaotization of solid medium structure in the course of melting 

is introduced. Thermodynamic and shear melting is described consistently. Dependences of friction 
force on temperature of lubricant, shear velocity of rubbing surfaces, and pressure upon surfaces are 
j^ analyzed. Within the framework of a simple tribological model the stick-slip mode of friction, when 

the lubricant periodically melts and solidifies, is described. The obtained results are qualitatively 



C/2 
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'^ compared with the experimental data. 
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^ 1 Introduction 

(N 

t^^ The boundary friction mode arising with lubricant thickness less than 10 atomic layers is widely in- 

O 

O vestigated so far ill. Experiments show that a thin layer of lubricant shows anomalous properties in 

Kj^ comparison with volume lubricants [2]. In particular, the intermittent motion {stick — slip) inherent 

S^ in a dry friction is observed f2|!3|. Such mode is explained as the solidification caused by compression 

c3 of rubbing surfaces followed by jump-like melting at a shear stress larger than the yield limit ("shear 

melting"). 

There are some phenomenological models, for example, thermodynamic [i], mechanistic [5] and 

synergetic Ig] ones, allowing to explain the experimental results partially. They are either deterministic 

J5}[6] or stochastic [7||8] in nature. Studies by using molecular dynamics methods are also known [9 12 



It turns out that the lubricant provides several kinetic modes with transitions leading to the stick — slip 



friction l2j. Such transitions are the first-order phase transitions 13 however between the states which 



are not true thermodynamic phases but the kinetic modes of friction. In work [T] three modes are found: 
shding at low shear velocities, regular stick — slip mode, and sliding at the high shear velocities. Their 



existence has been verified in numerous experiments Ufp, 14 . 

In work [g], within the framework of Lorenz model, to approximate viscoelastic medium, an approach 
has been developed according to which the transition of ultrathin lubricant film from solid- into liquid- 
like state results from the thermodynamic and shear melting. The processes due to self-organization 
of shear stress and strain fields and of lubricant temperature are described in view of additive noises 



of the specified quantities 15-17 and correlated temperature fluctuations '18'. The reasons for jump- 



like melting and hysteresis, observed experimentally 19-21 , are considered in works 22,23 . Here 



realization conditions of these features are determined with the deformation defect of the shear modulus 



taken into account. In work 24 , within the framework of the mentioned model the periodic stick — slip 
mode of friction is described which, however, has a stochastic component and can be realized only in 
the presence of fluctuations. The drawback of the model is that it does not consider the load applied 
to the surfaces of friction and that basic equations are derived by several approximations in [6]. 

In the present work, the theory based on the expansion of free energy of system into a power series 
of parameter / which, being squared, represents the excess volume arising from the formation of defect 
structure in the lubricant during melting. The liquid-like state is usually interpreted as a section of 
plastic flow in the loading diagram and characterized by the presence of defects in the lubricant |4]. 



Here, the approach based on the Landau theory of phase transitions is used 25 - 29 to describe strongly 
non-equilibrium processes, occurring at sliding of two rubbing solid surfaces separated by a layer of 
lubricant. 

2 Basic equations 

The melting of lubricants whose thickness is less than 10 molecular layers is judged by the increase 
in their volume [9] and diffusion factor [9|[l0|[30|[3l]. Since, of the two, the experimentally observable 
quantity is the volume, to describe the lubricant state we introduce parameter /, which, when squared, 
stands for the excess volume arising due to chaotization of structure of a solid medium under of melting. 
With the increase in / magnitude the density of defects in the lubricant is growing, and owing to their 
transport under the applied stresses the lubricant passes to a kinetic mode of plastic flow (liquid-like 
phase). 

Let's write down the dependence of the density of free energy $ on excess volume in the form of 
expansion in terms of parameter / 

$ = ^, + ^{Vff + ipof-\^lf + lv2f, (1) 



where $0) c, ipo, (pi, ^2 are expansion constants. 

We consider the dependences of elastic strains ef and lubricant temperature T on invariants for the 
smaller expansion powers only 



$0 = $S + 2M4r + /"(4) ' (2) 

ypo = ^l-\\{e%f -Ji{e\^^ -aT. (3) 

The first invariant represents the trace of strain tensor e^j = e\^ e\^ e\ and the second one is defined 



by the expression 32 



(4)' = (4)' - 2/2 = {e\ + el + elf - 1{e\e\ + e\el + e\e%) = {elf + {elf + (5^) 

According to ([!]), in the lubricant the elastic stresses appear 
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Equation (pi) can take the form of the effective Hooke law 



(^tj - '^Mf^lj + ^effStA 



ii"tj 



with the effective elastic parameterq^ 



which decrease with melting for increasing /. 
It is easy to show that 
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where n, r are normal and shear components of the stresses acting on lubricant from the side of rubbing 
surfaces^ Relationships pi), (10) represent the relationship between components of tensors and their 



invariants of the linear theory of elasticity 32 



Let's write down the kinetic equation for non-equilibrium parameter / in the form of the Landau- 
Khalatnikov equation 



^// 



a/' 



(11) 



^For p > [ijii, we put /Xg// = and for /^ > A/A, we take Ag// — 0. 



^Shear stress r is defined from expression (|6| for i 7^ j, i.e. 8ij = 0. When /ie// — the term t/ ^^ff hi (10), should 
be substituted by 2efj-, according to ([6|. 



where relaxation time Tf is introduced. In the exphcit form it is written as 

df ^ ^ . 2n2(A + /i)/ 



dt 



2' 



(12) 



with the last term showing the dependence of te| and (10) on /. 



We derive the equation connecting the relative velocity of shear for rubbing surfaces Vij and the 
elastic strains arising in the lubricant ef-. For this purpose we use the Debye approximation relating 
the elastic strain and the plastic one ef^- 4 : 



■pi ^ '^ 



(13) 



where r^ is the Maxwell relaxation time of internal stress. The total strain in a layer is defined as 



1-3 ^ U 



This strain fixes the velocity of motion of the top block according to equation 33 

he 
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(14) 



(15) 



where h is the thickness of lubricant film. The last three relationships give the expression for elastic 
components of shear strain: 

(16) 



Te^ij 



-£ • + 
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Further, for simplicity, a homogeneous system is considered, and in Eqs. ([I]), ( |I2| ) we put V = 0. 



3 Friction force 



In view of definitions (|3]), ([6]) - (10) the set of kinetic equations (12), (16) is closed and can be used for 
studying the kinetics of lubricant melting. In this section we consider the stationary modes of friction 



resulting from system evolution. According to (16), in the course of time, the stationary value of elastic 
component of shear strain is set in 

4o = 1^- (17) 

For finding the stationary states of all quantities it is necessary to solve numerically the evolution 
equation (12), using ([3]), ([6]) - (10) and determining the value of strain from (17). 

In experimental studies, the dependences of friction force on velocity of shear, thickness of lubricant 
layers, and load on the rubbing surfaces are often determined |l|j3| [I9}j2l] . Let's analyze how the 
lubricant temperature and velocity of shear affect the friction force. 

In the course of friction, both elastic (J% and viscous crji*^'^ stresses arise in the lubricant. The total 
stress in a layer is the sum of these two components 



(Ti 






The total friction force is defined in a standard way: 



F 



(JijA, 



(19) 



where A is the area of contacting surfaces. Viscous stresses in a layer are given by the formula 33 

vise Veff Vij 
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(20) 



where ?7e/-/ is the effective viscosity measured only experimentally, and in the boundary mode 



Veff ~ [Si 



IJJ ; 



(21) 



and for the majority of systems 7 = 2/3. In view of (15), (21) the expression for viscous stresses (20) 
is written in the form: 
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(22) 



After substitution of (18) and (22) in (19) we have the required expression for friction forcaj 
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(23) 



where af^ is defined by the formula m^ for z 7^ j. 

In experiments, the friction surfaces are often the atomically-smooth surfaces of mica, and lubricants 
are the quasispherical molecules of octamethylcyclotetrasiloxane (OMCTS) and linear chain molecules 
of tetradecane or hexadecane (2, 21 . The experimental conditions are as follows: the thickness of 
lubricant h ~ 10~^ m, the contact area A ~ 3 • 10~^ m^, load on the top surface of friction L = 
(2 -=- 60) ■ 10~^ N that corresponds to normal stresses n = —L/A = —(6.67 -^ 200) ■ 10^ Pa. The friction 
force F ~ (2 ^ 40) ■ 10~^ N. In the mentioned experiments it has been found that the lubricant melts at 
a temperature above the critical value T > Tco ~ 300 K, or at a velocity of shear V > Vc ^ 400 nm/s. 
These values may much differ depending on lubricant and geometry of the experiment. 

According to experimental data, the following values of the constants of the theory are chosen for 
the considered model: $* = 20 J/m^, A = 2 ■ 10" Pa, A = 10^ Pa, /i = 4.1 ■ 10" Pa, /i = 4 ■ 10" Pa, (^* = 
185 J/m3,(/?i = 570 J/m^,ip2 = 3200 J/m^a = 0.6 J-K'^/m^,h = 10"^ m, r/ = 1 Pa-s, r^ = 10"^ s. 
Note that the relaxation time of excess volume r/ has the dimension of viscosity. In fact, this means 
that the time of reaching the stationary friction mode increases with the growth of the effective viscosity 
of the lubricant. 

Figure [l^ illustrates the fact that the friction force decreases with lubricant temperature increase. 
Let's analyze the curve 1 in detail. At first, with temperature increase to point of phase transition 



^Here the sign function sgn(x) and absolute value of shear velocity \Vij\ are introduced, since Vij can accept negative 



values too. 
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Figure 1: Dependence of the stationary value of the total friction force Fij (niN) (23) on the temperature 



-7-10^ Pa: a) 



of lubricant T (K) and the velocity of shear Vij (nm/s) for 7 = 2/3, A = 3 ■ 10^^ m^, n ■ 

curves 1-4 correspond to the constant values of velocities of shear V^j=500 nni/s, 800 nni/s, 1100 nni/s, 

1550 nm/s; b) curves 1-3 correspond to the fixed values of temperature T=200 K, 255 K, 310 K. 

at T < Tco the excess volume is equal to zero, therefore the effective shear modulus 2/ie// ([T]) takes 



constant value. Since the shear velocity is invariable, in accordance with (17), the constant value of 



strain is realized, leading to constant shear stress ([6]). Thus, the both components in (23) do not change 
with temperature increase to a point of phase transition, and the friction force is constant. When the 
temperature becomes higher than the critical value T > Tco, the value of excess volume p increases 
abruptly, and the lubricant melts that leads to a sharp decrease in the total friction force. 

If after melting the temperature is increased, the excess volume increases resulting in reduction of the 
shear modulus and, as a consequence, in a smaller value of the total friction force. With temperature fall 
the lubricant solidifies when T = T^. Thus, the dependence is of hysteretic character that corresponds 
to the first-order phase transitions. With the increase in velocity of shear the lubricant melts at a smaller 
value of temperature. When the velocity exceeds the critical value, the lubricant is always liquid-like 
independent of temperature (curve 4), and the friction force decreases with temperature increase owing 
to the decrease in the shear modulus (the lubricant becomes more liquid). 

Thus, at small temperatures (T < T°) there is one zero minimum of the potential $(/) (solid- 
like lubricant). In the range of temperatures T^ < T < Tco, together with the specified minimum, two 
symmetric nonzero minima of $(/) coexistrl however, the system can not pass to the state corresponding 
to those minima since they are separated by maxima with the zero minimum. At further increase in 
temperature T > Tco the separating maxima disappear, and according to the mechanism of the first- 



*Sincc /^ is the observable quantity, these minima are equivalent in relation to friction mode. 



order phase transition the lubricant passes to the state corresponding to the minimum of potential with 
/ 7^ 0, i.e. it melts. If now the temperature is decreased, there appears a zero minimum and again 
the system can't pass to the respective state because of the presence of separating maxima. With their 
disappearance at T = T^ the lubricant solidifies abruptly. 

Figure ^p shows a somewhat different behavior. Here, in accordance with (23), at small shear 
velocities the lubricant is solid-like, and since / = the value of shear modulus is maximal. In such a 



mode an increase in the velocity increases the friction force (23) due to the increasing viscous stresses. 



At V^ > Vco the lubricant melts, and elastic shear stresses ^ decrease sharply and the total friction 
force decreases in a jump-like manner. With further increase in the velocity, the Fij increases owing 
to viscous component of the friction force, which is always increasing with the velocity. According 



to curve 3, in the liquid-like state the friction force (23) increases due to the increase in the velocity. 



With temperature increase the lubricant melts at lower shear velocities. Note, that the results of Fig. [TJd 



qualitatively coincide with the new map of friction for a boundary mode proposed in work 33 as a result 



of generalization of experimental data. At present, the dependences of friction force on temperature 
similar to those presented in fig. la are not measured experimentally. 



4 Stick-slip mode 

The dependences of figure [T] have been obtained for the fixed shear velocity of the top rubbing surface. 
However, the dynamic characteristics of tribological systems are defined by both the friction force, 
shown in the specified figure, and the properties of system as a whole. In particular, according to 
experiments, in the hysteresis region of the dependence in fig. [1} the interrupted mode of friction 
(stick — slip) (2| [SJ [7[]9| [2l| [3H [33] can be realized. The present work is devoted to definition of its 
features. The typical scheme of tribological system is presented in fig. |2] Here, the spring of rigidity k 
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Figure 2: The scheme of tribological system. 



is connected to the block of mass M to which additional normal load L is applied. The block is located 
on a smooth surface from which it is separated by a layer of lubricant of thickness h. The free end of 



the spring is brought in motion with a constant velocity Vq. Block motion initiates friction force F (23) 



that resists its displacement. Generally, for ultrathin layers of lubricants in a mode of boundary friction 



the velocities of block V and spring Vq do not coincide because of the oscillating character of force F 
leading to interrupted motion of the block. Such mode resembles dry friction without a lubricant. 
The equation of motion of the top block looks like l2lH|pp| 



MX = k {Vot -X)-F. (24) 



To calculate the time evolution of the friction force, the last equation is to be solved together with (12), 



(16), and the friction force defined from (23). However, the relaxation time of strain r^, may be put 



small in comparison with that of excess volume r/ owing to lubricant's thinness. Therefore, within the 



limits of approximation tj ^ r^, the two equations (24), (12) are solved jointly and the strain is defined 



from (17). 



The result of solution of the specified equations is shown in fig. |3| According to the figure, at first, 
the friction force monotonously increases, since the lubricant is solid-like, and shear velocity V increases. 
When the velocity exceeds critical value Vco, the lubricant melts and the friction force decreases, the 
shear velocity of the rubbing block increases, and it quickly moves to a large distance. As a result, the 
tension of spring and the velocity of shear decrease. When the velocity becomes less than that for the 
lubricant to be maintained in liquid-like state, it solidifies, and the friction force starts again increasing. 
The described process is periodically recurrent. It is worth noting that the velocity, with which the 
lubricant solidifies, does not coincide with the similar velocity shown in fig. [T| This is because of the 
sharp increase in the velocity of shear V at melting and corresponding increase in excess volume /^. 
According to ([T]), in this case the shear modulus becomes less than zero, but it need be considered zero 
and this deforms the potential form ([I]). Thus, in the presence of elastic strain (17), according to (|6|, 



the elastic stresses in the lubricant are equal to zero. As a result, the friction force decreases, and the 
lubricant flows. 

In fig. E^ there are dependences of the total friction force F (23), excess volume /^, and elastic 



component of the shear stress a?- I\6h on time under velocity Vq increase. Initially, the movement of 
the top sheared block (Vq = Vqi) leads to the increase in friction force at zero /^. When the elastic 
shear stresses reach the critical value, the shear melting of lubricant is initiated by the first-order phase 
transition mechanism. Now, parameter /^ increases abruptly, and elastic stresses become equal to 
zero. And again, there follows lubricant solidification since the relative shear velocity of the friction 
surfaces decreases (see fig. 3). When it solidifies completely, the elastic stresses appear there. And 
again, subsequent increase in the stresses leads to the increase in parameter /^ until it reaches the 
critical value needed for melting, and the process is repeated. As a result, the periodic interrupted 
{stick — slip) mode of melting/solidification becomes steady. With velocity increase to Vq = V02 the 
frequency of stick — slip peaks increases because at this velocity the critical value of stresses is reached 



^Further, for convenience we omit lower tensor designations since shear in one direction is considered. 




20 30 40 50 



70 t 



Figure 3: Dependence of the friction force F (mN), the shear velocity of rubbing surface V (nm/s), and 
its coordinate X (/im) on the time t (s) for n = -7 ■ 10^ Pa, M = 0.4 kg, k = 480 N/m, T = 250 K, 
Vo = 2500 nm/s. 



more rapidly. Accordingly, the lubricant melts more rapidly and in the same interval of time the system 
succeeds in making a greater quantity of melting/solidification transitions. With a higher increase in 
velocity Vq = Vqs, the frequency of stick — slip peaks decreases again due to the appearance of long 
kinetic sections F = const in the F{t) dependence. Note that in this mode, at melting, the excess 
volume /^ initially sharply increases owing to rapid increase in shear velocity of the top rubbing block 
V. A smaller value of excess volume /^ reached after abrupt initial slip of the rubbing surface at the 
expense of released mechanical potential energy of compressed spring corresponds to the stationary 
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Figure 4: Dependence of the friction force F (mN), the excess volume /^, and the elastic components 
of shear stress af^ (MPa) on the time t (s) for parameters of fig. pi normal stress n = —30 • 10^ Pa, and 
velocities of shear Vbi = 1800 nm/s, Vb2 = 3500 nm/s, Vqs = 4286 nm/s, Vo4 = 4290 nm/s. 



kinetic section. With further increase in velocity Vq = V04, the interrupted mode disappears, and the 
kinetic mode of liquid-like lubricant friction is reached, it is characterized by nonzero value of excess 
volume p and elastic stresses afj = 0. Thus, with the increase in velocity, the frequency of stick — slip 
peaks first increases, then it decreases due to the appearance of long kinetic sections, and when value 
of velocity exceeds Vq the stick — slip mode disappears. The described behavior is in a good agreement 
with experimental data [2]. 



10 



The influence of the external pressure, applied perpendicularly to the friction surfaces, on character 
of lubricant melting is often studied experimentally (2,31 . Such experiments show that pressure effect 
on the parameters of tribological systems is not trivial. For example, for a lubricant consisting of chain 
molecules of hexadecane the critical velocity of shear decreases with pressure growth, and for spherical 
molecules of OMCTS it increases, on the contrary |2|. Pressure influences the frequency and amplitude 
oi stick — slip transitions too [2]. Within the limits of our model, according to equation (12), the increase 
in load on the friction surface reduces the excess volume that should promote lubricant solidification. 

In figure [5} the time dependence of friction force is shown at various values of pressure acting 
to compress friction surfaces. At a temperature below the critical value (the top panel of figure) the 
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Figure 5: Dependence of the friction force F (mN) on the time t (s) for parameters of fig. pi shear velocity 
Vq = 1600 nm/s, and external normal load rii = —7 ■ 10^ Pa, n2 = —70 ■ 10^ Pa, n^ = —100 ■ 10^ Pa, 
n4 = —180 ■ 10^ Pa. The top panel corresponds to the temperature T = 235 K, the bottom one to 
T = 540 K. 

interrupted mode of friction is realized, and with pressure growth the amplitude of stick — slip transitions 
and values of kinetic and static friction force increase, as to the frequency of transitions, it decreases. 
At the pressure, corresponding to normal stress n = n^, the stick — slip mode is not realized. However, 
the kinetic mode, corresponding to lubricant melting, does not set in and lubricant solidifies because 
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of walls squeezing. As a result, the lubricant can't melt, the friction force F is high, it corresponds to 
solid-like lubricant and zero value of excess volume /^, as the squeezing of walls promotes the formation 
of long-range order of atoms in the lubricant. In the bottom figure panel, the dependence is for the 
elevated temperature T. Here we have the kinetic mode of friction corresponding to a small value of the 
friction force and to nonzero excess volume p . With pressure growth the stick — slip mode is expected, 
and with a higher loading the full solidification of the lubricant is to be observed the same as in the 
top panel of the figure for n = n^. Thus, three modes of friction have been revealed: 1) the kinetic 
mode in which the lubricant is always liquid-like; 2) the interrupted mode corresponding to periodic 
melting/solidification; 3) the mode of dry friction characterized by larger value of the friction force and 
solid-like lubricant structure. Similar modes are found also within the framework of stochastic model [T] . 
In figure [6| the dependence of friction force on lubricant temperature is shown. Temperature rise 
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Figure 6: Dependence of the friction force F (mN) on the time t (s) for parameters of fig. pi velocity 
Vq = 3850 nm, and lubricant temperatures Ti = 150 K, T2 = 300 K, T3 = 650 K, T4 = 800 K. 



leads to the decrease in the amplitude of friction force oscillations and the change of the frequency of 
phase transitions. For T = T4 the mode of sliding arises characterized by constant values of kinetic 
friction force and shear velocity of the rubbing block. Thus, the temperature growth promotes lubricant 
melting. Experiments with a similar study of the temperature effect are unknown to us and the given 
dependence is therefore predicting. 

5 Conclusion 

The proposed theory allows to describe the effects observed at the melting of ultrathin lubricant film in 
a boundary friction mode. The consistent consideration of thermodynamic and shear melting has been 
carried out. The dependences of friction force on shear velocity and temperature, as well as on pressure 
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upon surfaces are analyzed. At high lubricant temperature the shear melting is realized at a smaller 
value of shear velocity (shear stresses), and with a still higher increase in temperature the lubricant 
melts even at zero velocity of shear. The model takes the effect of temperature, shear melting, and 
pressure into consideration. These are the major factors investigated experimentally. 

Within the framework of the proposed theory a simple tribological system has been studied, and 
the time dependences of friction force are obtained for the increasing shear velocity, pressure, and 
temperature. It is shown that in the system, in a wide range of parameters, the experimentally observable 
intermittent mode of friction is realized. The obtained results qualitatively coincide with the known 
experimental data. Since the model is quantitative, it may be modified to describe concrete experiments. 

The work has been done under financial support of the Fundamental Researches State Fund of 
Ukraine (Grant $28/443-2009). 
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